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We use the numerical renormalization group method to 
study the 0(3)-symmetric version of the impurity Anderson 
model of Coleman and Schofield. This model is of general 
interest because it displays both Fermi liquid and non-Fermi 
liquid behaviour, and in the large U limit can be related to the 
compactified two channel Kondo model of Coleman, Ioffe and 
Tsvelik. We calculate the thermodynamics for a parameter 
range which covers the full range of behaviour of the model. 
We find a non-Fermi liquid fixed point in the isotropic case 
which is unstable with respect to channel anisotropy. 
PACS Number: 7520H 



I. INTRODUCTION 

The unusual experimental results observed in many as- 
pects of the behaviour of the high temperature supercon- 
ductors and some heavy fermion systems Q, has led to 
the conjecture that this behaviour can best be explained 
by some form of non-Fermi liquid theory. For strong cor- 
relation models that have been put forward to describe 
these systems it is very difficult to test the stability of 
the Fermi liquid state due to the strong local interac- 
tion U. Perturbative methods are unreliable for large 
values of U, and non-perturbative techniques are in gen- 
eral difficult to apply and subject to uncertainties in the 
approximations that have to be used. With direct numer- 
ical methods such as Monte Carlo calculations there are 
very severe size limitations (and also the 'sign problem') 
that make it difficult to obtain results for excitations for 
the thermal energy scale that are required to examine 
these issues. It is of some interest, therefore, to study 
strong correlation impurity models that display similar 
types of non-Fermi liquid behaviour as accurate methods 
for predicting the behaviour of these systems have been 
developed, including many exact solutions. 

The non-Fermi liquid behaviour of the impurity two 
channel Kondo model has already been studied exten- 
sively and has been used as a basis for explaining the 
behaviour of certain heavy fermion compounds It 
has also been used to interpret the experimental results 
for certain two level tunnelling systems || . 

Recently a modified form of the Anderson impurity 
model has been formulated by Coleman and Schofield ||] 
which also displays non- Fermi liquid behaviour. In the 
large U limit the model can be related to a form of the 
two channel Kondo model Q, and for small U there is a 
parameter regime where it displays marginal Fermi liquid 
behaviour ||||. This model is simple enough so that its 



behaviour can be studied using the numerical renormal- 
ization group techniques (NRG) that were used to ob- 
tain definitive results for the standard Anderson model 
0. This modified Anderson model, the 0(3)-symmetric 
Anderson model (0(3)-AM), has the form, 

#0(3)-AM = £>/J/a + U f\M\h 
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Apart from the term with the anomalous hybridization 
matrix element V a this is the standard Anderson model 
for a localized f level hybridized with a conduction band 
and an on-site interaction U. The model is taken to be 
particle-hole symmetric with £f = —U/2 and £/- = — e_fc. 
Due to the antisymmetry of the dispersion in -ffo(3)-AMi 
the hybridization V(e) and the anomalous hybridization 
V a (e) have to be antisymmetric (and are here chosen to 
be constant). 



V(s) = 



V : e > 
-V : e < 



14(e) = 



F a : £>0 
-K : e<0 



(2) 



The extra anomalous hybridization term is one in which 
neither particles nor spin is conserved. It reduces the 
symmetry of the standard Anderson model at particle- 
hole symmetry from 0(4) to 0(3) (the 0(3)-symmetry is 
best seen in the Majorana Fermion formulation ||||]). 

The motivation for its introduction is that in the large 
U limit the model can be transformed by a Schrieffer- 
Wolff transformation Q into a localized s-d type of model 
of the form 
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£kc\ a c ka + S ■ [J x s(0) + J 2 f(0)] . (3) 



The operators s(0) are the spin operators for the conduc- 
tion electron states defined in terms of the conduction 
electron creation and annihilation operators (c\ a , c k a) in 
the usual way, 
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The r(0) operators are isospin operators for the con- 
duction electrons and these are defined by 
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(5) 



Both sets of operators satisfy the usual SU(2) commu- 
tation relations. The spin and isospin of the conduc- 
tion electrons are coupled to the impurity spin S indi- 
vidually with coupling constants J\ and Ji , respectively. 
From the Schrieffer- Wolff transformation Q these cou- 
pling constants are given by 



j _ W{V - V a ) T _ WV a 
Ji- jj , -h- — 



(6) 



This localized model eq. (3), which was introduced by 
Coleman, Ioffe and Tsvelik ]^|, has a form very similar 
to the two channel Kondo model and it has been referred 
to as the compactified two channel Kondo model or (r- 
cr)-model. Though it differs from the two channel Kondo 
model in that the impurity spin cannot be overscreened 
by the conduction electrons (the states with spin and 
the states with isopin are mutually exclusive and so the 
cannot both be used to screen the impurity at the same 
time) it has been argued that it has the same low energy 
fixed point as the two channel Kondo model ||. The 
channel isotropic case (Ji = J2 in eq. (6)), where we 
expect non-Fermi liquid behaviour, corresponds to V a = 
V/2. 

It is clear that the general 0(3)- AM has a rich range of 
behaviour and may provide theoretical insights into the 
controversial questions as the non- Fermi liquid behaviour 
observed in some strongly correlated systems. Due to the 
particle non-conserving terms it is unlikely that the Bethe 
ansatz technique, which has been used to solve many of 
the strong correlation impurity models, can provide a 
solution for this case. 

The numerical rcnormalization group method, how- 
ever, as originally developed by Wilson for the Kondo 
problem and further used by Krishnamurthy et al. 
jjj to obtain definite results for the Anderson model in 
all parameter regimes can be used. In this paper we apply 
this approach and give details of the calculation and re- 
sults for the thermodynamic behaviour (a brief summary 
of this work is given in [^)). The only slight complica- 
tion in applying the NRG to the 0(3)- AM is that the 
many-body states obtained in the iterative diagonaliza- 
tion cannot be classified in terms of quantum numbers 
associated with the particle number and spin individu- 
ally. There is a total spin operator T, which is the sum 
of the spin and isospin, which does commute with the 
Hamiltonian. It is defined by 



where the sum runs over the impurity and conduction 
electron states. For perturbational calculations it is use- 
ful to express the Hamiltonian in terms of Majorana 
fermion operators In the Majorana fermion form 

it is clear that the non-Fermi liquid situations which de- 
velop for V a = V/2 are due to the fact that one of the 
Majorana fermions associated with the impurity is unhy- 
bridizcd. This gives a local zero mode and the scatter- 
ing with the conduction electrons via the U term leads 
to singularities and the breakdown of Fermi liquid the- 
ory. There is no advantage, however, in the numerical 
renormalization group approach in using the Majorana 
fermion representation so we retain the model in the form 
eq. (|l|). We do need, however, to re-express the Hamil- 
tonian along standard lines for the numerical renormal- 
ization group approach and introduce a discrete basis for 
the conduction electrons for the iterative sequence of di- 
agonalizations that is used. 



II. DETAILS OF THE NUMERICAL 
RENORMALIZATION GROUP ANALYSIS 

In order to apply the NRG method on the 0(3)-AM 
eq. ([!]), its Hamiltonian has to be cast in a semi- infinite 
chain form. This mapping is described in detail in (Jt|. 
The logarithmic discretization (with the discretization 
parameter A > 1), the definition of the discrete opera- 
tors a npCT and b np(T and the transformation of the free 
conduction electron term are identical to the treatment 
of the standard single impurity Anderson model (SIAM) 
0. The p ^ fermions again do not couple directly to 
the impurity and are neglected. 

Due to the antisymmetry of V{e) 1 there is a minus sign 
in the definition of the conduction electron operator at 
the impurity site 
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In contrast to the combination a na + b nr7 one obtains in 
the standard case. 

The discretized version of -ffo(3)-AM then reads 

#0(3)-AM = 5>f/^ + UfiMlh 

+ Vj2(ftgo, + glu) (9) 
- K (glji + 501/1 + /Jsoj + fUli] 
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The minus sign in (8) does not affect the semi-infinite 
(7) chain form of the conduction electron part 
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with 



*n =2(l+A- 1 )(l-A" n - 1 ) 

»-2n-;IA _1 / 2 / 1 _ A-2n-3\ — V 2 A-n/2 



x (l-A-^-^^l-A- 2 "- 3 ) 



(11) 



As in the standard case we diagonalize iteratively a 
sequence of Hamiltonians Hjq with 
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(13) 



so that 



H, 



0(3)-AM 



= lim A-^-^^H 



N- 



N 



(14) 



At this point it is convenient to generalize the Hamilto- 
nian (9) by an anomalous hybridization between all sites 
of the semi-infinite chain. Together with the V a -term we 
have 



Ha — ^ H a ri , 



(15) 



9n+ll9nl 



\n+l 



+ g^Qn+n + (-i) n+1 9ld + n 



Qn+Ugnl 



(16) 



with i aj _i = V a . The (— l) n+1 factors are due to the 
requirement that T 2 commutes with H. Obviously, all 
t a . n are zero for n > 0, but the advantage of this gen- 
eralization is that we can start the iterative procedure 
by adding site to the impurity. Otherwise the two-site 
problem would have to be the starting point. 

Both particle number Q and spin S are no more 
conserved quantities of the Hamiltonian (9) as soon as 
V a =/= 0. The only conserved quantity is the total 'spin 
plus isospin' 



T = S + t 



with 



T = T[ 



71=0 
CO 

s = sf + s,, 



(18) 
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and 



t+ - f-ir+V o ] 

T n — \ L ) 9n]9n[ > 

Tn = (-l)" +1 3«i5nT 

Tnz = g (fl£ T 5nt + 9ii9nl ~ 1 



(20) 



The n-dependent definition of r„ is necessary to fulfill 
the condition [T + /~,H] = 0. 

The starting point for the iterative diagonalization is 
the uncoupled impurity (labeled by —1). Writing the 
states in the form \T,T Z , r) we get 
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Supposed that, at the iVth step of the iteration proce- 
dure, the Schrodinger equation 

H N \T, T z , r) N = E N (T,r)\T,T z ,r) N (22) 

has been solved, we can set up the Hamiltonian matrix 
for the (N + l)st step. First we define a basis for the 
enlarged system 

\T,T z ,r) N , 



\T,T z ,r;0) 

\T,T z ,r;]} = g ] N+1] \T,T Zl r) N , 

\T,T z ,r;i) = g i N+n \T,T z ,r) N , 

\T,T z ,r;U) = g ] N+1 ^ N+n \T,T Zl r) 



(23) 
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These states can be combined to give eigenstates of T 2 
and T 7 



\T,T z ,r;l) 



N+l 



{-rfJW\T~hT z ~\,r-U) 
1 '• 'T-±,T z + i,r;0) , 



2T 



\T, T Z: r; 2) N+1 



2T 
T-Tz 



T+T z rp _ 1 rp _ l *\ 

-t 2' - Lz 2' ' ' 1 / 

r-ir i + i,r;i) , 



2T 



(17) |T,T 2 ,r;3) w+1 = 



2T+2 - t ^2' z 2' ' ' 1 / 

T+T. + l | T , 1 rp 1 ,\ 
2T+2 | J -^2' Z '2' '' 1 -/ 



1^, T z , r; 4) w+1 — 



T-T z + 1 
2T+2 



T+i,T 2 -|,r;Tl) 



<-V N yf 



N /T+7V+1 | T , 1 rp , 1 M 
2T+2 | J ' 2' z ' 2> ' ' U / 



n=0 



(24) 
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The Hamiltonian matrix 

H T (ri,r'j) = N+1 (T,T z ,r;i\ H N+1 \T,T z ,r';j) 



N+1 

(25) 



now consists of three parts: 

H T (ri,r'j) = VX n+1 (T,T z ,r;i\ H N \T,T z ,r';j) N ^ 
+ k N ' 2 Y,tN N+i(T,T z ,r;i\X\T,T z ,r';j) m 

a 

+ A N/2 



The matrix elements of eq. (26) are given in the appendix. 

We now have all the information to set up the ma- 
trix Hx(ri,r' j). The following steps are completely 
analogous to the standard case. By diagonalization of 
HT(ri,r'j) the Schrodinger equation of the (N + l)st 
step is solved. Only a limited number of the lowest lying 
states (500 in our calculations) are kept and used to set 
up the Hamiltonian matrix for the next step. 



+1 (T,T z ,r;i\H aN \T,T z ,r';j) N+1 . (26) 



N 



with X = glf^gN+ia + 9^ N+la 9Na- The first part simply 
gives the eigenstates of the previous iteration 

N+i (T,T z ,r;i\ H N \T,T z ,r';j) N+1 = 5 rr >5 i3 



x < 



f E N (T-l/2,r) 
E N (T- 1/2, r) 
E N {T + l/2,r) 
{ E N (T + l/2,r) 



i = 1 

i = 2 

i = 3 

i = 4 



(27) 



To make use of the Wigner-Eckart Theorem 

1 



{T,T z ,r-i\V*\T',T' z ,r';i) = 



(T',r'\\V k \\T,r) 



V2T' + 1 v ~ '' 11-9 
x(T,T z ,k,q\T'X) (28) 

for the calculation of the second and third part, we have 
to define irreducible tensor operators according to 



[T* V k ] _ = ^/k(k + l)-q(q±l) V q k ±1 , 
[T z ,V q k ]_ = qV q k . (29) 

In the standard case (where we have the spin S instead 
of T in eq. (29)) the tensor operators are just the anni- 
hilation operators f± and /f . Here we have 



V° - 



i((-ir<4+5ni) , 



n,0 ^ 

^o = ^((-l)" + V U +5«l) , 



(30) 



The and g na can be written in a form which only 
depends on the operators V„ and V® 



5i T 



^=(-l)" +1 [T+.^o-^o]. 
2li = ^(-l) n+ MKV<o) » 



(31) 



III. DISCUSSION OF THE FIXED POINTS 

The various possible fixed points of the 0(3)-AM can 
be identified by investigating the flow diagrams for the 
lowest lying eigenenergies. Fig. 1 shows a typical flow 
diagram for V a = (which corresponds to the standard 
SIAM). The solid (dashed) lines belong to states with 
T = 0.5 (T = 1.5). The ground state has T = 0.5. This 
flow diagram is similar to Fig. 6 in with the difference 
that in § the energies are labeled by particle number 
Q and spin 5*. The system flows from the Free Orbital 
to the Strong Coupling fixed point. The influence of the 
Local Moment fixed point is small for the parameters 
used in Fig. 1 but increases for larger values of U. 

In Fig. 2 we present a flow diagram for the isotropic 
case 14 = V/2. We clearly find a new fixed point, which 
we call the non-Fermi liquid fixed point (see the discus- 
sion in the following sections). The non- Fermi liquid 
fixed point is unstable for small deviations from isotropy 
as can be seen in Fig. 3 for V a = 0.498V. Up to iter- 
ation number 30 the flow diagrams of Fig. 2 and 3 are 
almost identical but after a crossover regime (iterations 
30 to 50) the system approaches the Fermi liquid fixed 
point of the single channel case (see Fig. 1). The inter- 
pretation of this as a Fermi liquid fixed point can also 
be found in Fabrizio et al. |lC|| , who investigated the two 
channel model by a bosonization method. 

The static properties are calculated at each iteration 
step N for the temperature 



T N = Ia-^ 1 )/ 2 





(32) 



For (3 we take the value (3 = 0.46 (see the discussion in 
) . Due to this connection of the number of iterations 
and the temperature we can see that for V& — 0.498V 
there exists a high temperature regime in which the sys- 
tem behaves similarly to the isotropic case and a low tem- 
perature regime with a Fermi liquid ground state. This 
does not automatically mean that we have non- Fermi liq- 
uid properties in the high temperature regime because 
the non-Fermi liquid physics itself dominates only in the 
very low temperature regime of the isotropic case. 
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IV. RESULTS FOR ENTROPY AND SPECIFIC 
HEAT 

At each step of the iterative diagonalization the impu- 
rity contribution to the Free Energy is calculated by 



F{T N ) 



-k B T N 



fiE" 



(33) 



with Tjv given by equation (|32|). The Ei and E® are the 
scaled energies of the full system and the system without 
impurity, respectively. The entropy S(T) and the spe- 
cific heat coefficient 7(T) = C(T)/T are calculated by 
numerically differentiating the Free Energy 



S(T) = ~F(T) , 



dT 2 



F{T) 



(34) 
(35) 



It is useful to take the values of F(Tpf) only at odd (or 
even) N. One has to be careful about the groundstate 
energy -Eg.jv which is substracted at each iteration step 
and has to be added again to give the correct derivative. 
In the formula for the entropy 



S[l(T N + T N - 2 )\ 



F{T N ) + AE- F{T N _ 2 ) 

Tn — IjV-2 



the correction term AE takes the form 
AE = A"^- 1 )/ 2 f VAE.jf-! + E 



(36) 



(37) 



It consists of two terms because we compare the free en- 
ergy between two steps of the iteration procedure. 
For the specific heat coefficient we use 



7 



S[$(T N +Z 



N-2 



\T N _ 2 + i(T N 
)]-S [±(T W _ S 



T N -i)} = 

f Tjv-4)] 



(38) 



Results for the entropy are shown in Fig. 4. For V a = 
the entropy vanishes for T — > 0. In the isotropic case 
V a = V/2 we find a zero point entropy of 1/2 In 2 (within 
the numerical accuracy). For any < V a < V/2 we 
find the behaviour that has already been discussed in the 
previous section. The entropy is identical to the result 
of the isotropic case in a high temperature regime, but it 
tends to zero for T — ► 0. 

The specific heat coefficient j(T) is plotted in Fig. 5 
for V a = 0, V a = V/2 and V a = 0.43V. In the single 
band limit we find the usual Fermi liquid behaviour with 
a constant 7 at low temperatures, while in the isotropic 
case there is a very clear logarithmic dependence. For 
any V a V/2, j(T) approaches a constant value (7) for 
T — > 0. The energy scale T* (defined as I/7) depends 
quadratically on (V — 214) (not plotted here). This result 



corresponds to the calculations by Pang and Cox |Tl[] who 
found a T* oc (AJ) 2 dependence for the two channel 
Kondo model with channel anisotropy. 

In the isotropic case we define the energy scale T\ = 
1/c with c the prefactor of the InT term in 7. In the 
large U limit we find Tj oc exp(— const • U/irA) (see Fig. 
6, A = l/27rV 2 ) which corresponds to the expression 
for the Kondo temperature in the two channel Kondo 
model (Ik oc exp(— const/ J), J oc A/U). For U = the 
prefactor c vanishes and 7 approaches a constant value 
proportional to T _1 , with T a second energy scale not 
present in the large U limit. In the small [/-regime the 
energy scale 3] diverges as U~ 2 (see Fig. 6). This result 
is in agreement with second order perturbation theory 
around U = Ml. 



V. RESULTS FOR 'SPIN PLUS ISOSPIN' 
SUSCEPTIBILITY AND WILSON RATIO 

The most relevant susceptibility to calculate is \' = 
Xo-,imp + Xcr,c + Xt,c, which includes spin and isospin of 
the conduction electrons but only the spin of the impu- 
rity, x' is related to the spin susceptibility of the two 
channel Kondo model where x<r,c and Xr,c represent the 
spin susceptibility of channel one and two, respectively. 
The quantity x' ■> however, is not conserved and therefore 
difficult to calculate within the NRG (additional matrix 
elements have to be calculated iteratively) . 

On the other hand, the 'spin plus isospin' susceptibility 
Xt = X<r,imp + Xr,im P + X<x,c + Xr, c is conserved and hence 
can be calculated from the energy states only using the 
relation 



XT 



i T 2\ = T z exp 
^ z > Y,i exp(-/3S, 



<-0Ej) E^ z 2 exp(-/3£;°) 

-Be,) Y,i*M-M) 



(39) 



(40) 



with T z the ^-component of the total 'spin plus isospin' 
T. 

In the limit of large U the susceptibility Xr.imp vanishes 
because charge fluctuations (described by the isospin r) 
are frozen out at the impurity site. That means that, 
although we are studying the susceptibility \t for all 
values of U, it corresponds to \' ln the limit where 0(3)- 
AM and (T-er)-model are related via the Schrieffer- Wolff 
transformation. 

In the isotropic case (V a = V/2) we find xt oc InT 
while for any V a ^ V/2, \t approaches a constant value 
for T — > (not plotted here) . 

The Wilson ratio xt/i is plotted in Fig. 7 for both 
V a = and V a = V/2. Also shown is the result for 
Xs/l (xs + Xa.c) calculated with the standard 

program for the SIAM. In contrast to xs/l the 'spin 
plus isospin' Wilson ratio is independent of U (within 
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the numerical accuracy). For V a = V/2 we find a Wilson 
ratio which is consistent with the value of 8/3 of the two 
channel Kondo model, within the numerical accuracy. 

For < V a < V/2 the Wilson ratio is a nonuniversal 
function of V a (see Fig. 8). Due to the vanishing energy 
scale in this limit the numerical calculation of \t and 7 
gets more and more inaccurate. However, it is clear that 
R decreases rapidly as V a — > V/2 (note that the xt/i~ 
ratio is discontinuous at this point and for V a — V/2 the 
Wilson ratio takes the value 8/3). 



the same for the 0(3)-AM and the two channel Kondo 
model. The generalization of the NRG for the calculation 
of these quantities is in progress. 
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us (R.B.) was supported by a grant from the Deutsche 
Forschungsgemeinschaft, grant No. Bu965-1/1 and we 
also acknowledge support from the EPSRC (grant No. 
GR/J85349). 



VI. SUMMARY 



To summarize, we have discussed a numerical renor- 
malization group approach to the 0(3)-symmetric An- 
derson model that can be related to the (r-cr)-model 
introduced by Coleman et al. || via a Schrieffer- Wolff 
transformation. Both spin and charge are not conserved 
in the 0(3)-symmetric model so we have to use the 
only conserved quantity 'spin plus isospin' to classify the 
eigenstates in the iterative diagonalization procedure. 

The flow diagrams show a non-Fermi liquid fixed point 
for V a = V/2 (corresponding to J% — J2 in the (r-cr)- 
model) which is unstable for any small deviations from 
the isotropic case. All static properties calculated for 
y a = V/2 show the non- Fermi liquid behaviour of the 
two channel Kondo model: 

• a zero point entropy of i In 2, 

• a logarithmic temperature dependence of the spe- 
cific heat C(T) ocTlnT , 

• a logarithmic temperature dependence of the 'spin 
plus isospin' susceptibility Xt{T) oc InT , 

• a Wilson ratio of xt/7 = 8/3 independent of U. 

The relation of these results to those of the Conformal 
Field Theory and Perturbation Theory and the interpre- 
tation of the non-Fermi liquid fixed point in terms of 
Majorana Fcrmions will be dealt with in a separate pub- 
lication. 

For V a ^ V/2 the model shows the conventional Fermi 
liquid behaviour similar to the standard SIAM. There 
is no indication of a new non-Fermi liquid fixed point 
as conjectured in 0). We find an energy scale T* cx 
(Ji — J2) 2 which corresponds to the result of Pang and 
Cox 0|. This energy scale should correspond to T a , the 
anisotropic energy scale in |Q, but there T a oc (Ji — Ji). 
The difference is probably due to the different cut-off 
scheme used (giving a different prefactor to the exponen- 
tial). 

It is interesting to see whether the behaviour of 
non-conserved quantities (such as charge- and spin- 
susceptibility) , dynamic and transport properties are also 
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APPENDIX A: 



In this appendix we want to list the results for the matrix elements necessary to set up the Hamiltonian matrix 
(25). For the matrix element 

N+i(T,T z ,r;i\^ NtT g N+la +gl f+lrT g Na ^\T,T z ,r';j) N+1 (Al) 



we get the following results 
i = 2, j = 1: 



i _ 9T 1 

n(T - 1/2, r\\V^ Q \\T - 1/2, r') N + — = N (T 1/2, r||y# i0 ||T - 1/2, r') N ; (A2) 



4VT^/T 2 - 1/4" X -/-"ir^oir ' / iv 2\/r w- 
i = 2, j = 4: 

1 



i = 3, j = 1: 



* = 3, j = 4: 



( ~ 1)Jy ^+T Af ^ T " 1/2 ' r H^°H T + 1/2 ' r V (A3) 



^^^(r + 1/2, r|| V* |T - 1/2, r') N , (A4) 



i-ir^^^ N {T + i/2 A yhAT^y) N m 



+ (-l) N ^L^jN(T+l/2,r\\vl \\T+l/2,r>) N . (A6) 

The remaining matrix elements (i = 1, j = 2 etc.) follow from the hermiticity of H(ri,r'j). The anomalous part 
in the Hamiltonian matrix 



N+1 (T, T z , r; i\ H aN \T, T z , r';j) N+1 (A7) 



contains the contributions 
i = 2, j = 1: 



i = 3, j = 4: 



1 

7f 



(T-l/2,r|y° i0 ||T-l/2,r% ; (A8) 



(-1)^-^=^^+1/2^1^0117+1/2^% . (A9) 

The next step is to relate the reduced matrix elements n < l^jvoll >N an< ^ N < ll^/vol >n to the unitary matrices 
Ut of the previous iteration (the Ut diagonalize the Hamiltonian matrices Ht)- The results are 



(t,w\\v^ \\t + i,w') n = 



{ _ 1)N+ iVE+^^±l £ [Ur{Wt r4 )£/ T+1 {w >, r2) + U T (w, r3)U T+1 («,', rl)] , 



(A10) 



N (T,w\\V}j fi \\T, w') 



N 



7 



X ' 12 " lVl ^ l -^2[{-l) N U T {w,r2)U T {w' ,rl) -U T {w,rl)U T {w\r2)\ 



2T 

i 

+ ^^^^^[U T (w,r3)U T (w\rA) + (-l) N + 1 U T (w,r4)U T (w',r3)] , (All) 
jv(T,«;||l4 i0 ||T- l,w') N = 

V '- T 4 1 v : ' / 1 ^ [C/ T ( U ;,r2)C/ T _ 1 K,r4) + t^(t«,n)i;r-iK,r3)] , (A12) 



2T 



JV (T,t«||V^ i0 ||T,i« , ) JV 



T+ \J2 [(-l) N U T (w,r2)U T (w\rl) - U T (w, r3)U T (w' , r4) 

r 

+C/ T (w,rl)C/ T (w;',r2) + (-l) Ar+1 J7 T (w, r4)C/ T (w', r3)] (A13) 
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0.0 1 1 ==- ■ 1 • 1 

0.0 20.0 40.0 60.0 

N 

FIG. 1. Flow diagram for e f = -0.000714, U = -2e f , 
V = 0.01414 and V a = 0. Solid and dashed lines belong 
to T = 0.5 and T — 1.5, respectively. The system flows to the 
Fermi liquid fixed point. 



0.0 1 ' ■ ' • 1 

0.0 20.0 40.0 60.0 

N 

FIG. 3. Flow diagram for e f = -0.000714, U = -2e f , 
V = 0.01414 and V a = 0.00704. Solid and dashed lines be- 
long to T = 0.5 and T = 1.5, respectively. After a crossover 
regime, the system flows to the Fermi liquid fixed point. 




60.0 



N 

FIG. 2. Flow diagram for e f = -0.000714, U = -2e f , 
V = 0.01414 and K = V/2. Solid and dashed lines belong to 
T = 0.5 and T = 1.5, respectively. The system flows to the 
non-Fermi liquid fixed point. 



2.0 




FIG. 4. Entropy for e f = -0.000714, ?7 = -2e f , 
V = 0.01414 and V a = 0, 14 = V/2, and various values of 
< Ki < V/2 . In the isotropic case we find the residual 
entropy of 1/2 In 2. 
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FIG. 5. Specific heat coefficient j(T) for e f = -0.000714, 
U = -2e f , V = 0.01414 and V a = (solid line), V a = 0.0061 
(dotted line) and V a = V/2 (dashed line). 



-15.0 



-20.0 




-10.0 



2.0 3.0 
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FIG. 6. Energy scale in the isotropic case 14 = V/2. 
For small U, T ; is proportional to U~ 2 (insert) but goes as 
exp(— const • U/nA) for large U. 
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FIG. 7. Va-dependence of the Wilson Ratio R for 
et = -0.000714, U = -2e f and V = 0.01414. The V a ^0-case 
corresponds to the ordinary SIAM. The numerical calculation 
of R becomes more and more inaccurate as V a approaches 
V/2, (vertical line) but the data indicate that R vanishes in 
this limit. 
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FIG. 8. Wilson Ratio X t/i for V a = and V a = V/2 . Also 
shown is the Wilson Ratio xs/l for the standard SIAM. 
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